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基于 Canopy 聚 类 的 噪声 自 适应 模糊 C- 均 值 算法 
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摘 要 : 针对 局 部 空间 信息 的 模糊 C- 均 值 算法 (WFLICM) 中 空间 影响 因子 容易 受到 噪声 影响 出 现 错误 标识 的 问题 ， 
提出 一 种 融合 局 部 和 非 局 部 空间 信息 的 模糊 C- 均 值 聚 类 图 像 分 割 算 法 (NLWFLICM)， 在 WFLICM 算法 的 模糊 影响 因 
子 中 引入 非 局 部 空间 信息 ， 根 据 噪声 程度 自 适 应 地 设置 局 部 和 非 局 部 信息 权重 ， 并 重新 标记 中 心 点 的 模糊 影响 因子 。 
实验 结果 表明 , NLWFLICM 算法 具有 比 WFLICM 算法 更 强 的 鲁 棒 性 和 自 适 应 性 ,并 在 一 定 程度 上 提高 了 WFLICM 算 
法 对 含有 大 量 噪声 图 像 进 行 分 割 的 鲁 棒 性 , 同时 保留 了 图 像 的 纹理 。 为 了 提高 算法 的 聚 类 性 能 和 收敛 速度 , 结合 Canopy 
算法 能 够 快速 对 数据 进行 粗 聚 类 的 优点 , 提出 基于 Canopy 聚 类 与 非 局 部 空间 信息 的 FCM 图 像 分 割 改进 算法 (Canopy- 
NLWFLICM)， 可 以 在 NLWFLICM 算法 聚 类 前 ， 对 聚 类 中 心 进行 预 处 理 ， 从 而 提高 收敛 速度 和 图 像 分 割 精 度 。 
关键 词 : 聚 类 算法 ; Canopy 算法 ; 模糊 C- 均 值 算法 ; 局 部 和 非 局 部 空间 信息 
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Improved fuzzy C-means clustering algorithm based on canopy clustering 


Chen Kai, Chen Xiuhong, Sun Huiqiang 
(School of Digital Media, Jianenan University, Wuxi Jiangsu 214122, China) 


Abstract: Aiming at the problem that the spatial influence factors are easily misidentified by noise in the fuzzy C-means 
algorithm (WFLICM) for local spatial information, this paper proposed a fuzzy C-means clustering algorithm for image 
segmentation (NLWFLICM) based on local and non-local spatial information. It introduced the non-local spatial information 
into the fuzzy influencing factor of WFLICM algorithm, the weight of local and non-local information is adaptively set according 
to the noise level, and the fuzzy influence factors of the central point are re-marked. The experimental results show that the 
NLWFLICM algorithm is more robust and adaptive than the WFLICM algorithm, and improves the robustness of the WFLICM 
algorithm to a large extent, while preserving the image texture. In order to improve the clustering performance and convergence 
speed of the algorithm, combined with the advantages of Canopy algorithm for fast clustering of data, this paper proposes an 
improved algorithm for FCM image segmentation based on Canopy clustering and non-local spatial information (Canopy- 
NLWFLICM) before clustering algorithm. This can improve the convergence speed and image segmentation accuracy. 
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果 不 理想 。 为 此 ， 国 内 外 学 者 研究 了 许多 基于 局 部 空间 信息 的 


改进 FCM 算法 。Ahmed 等 人 外 在 FCM 的 目标 函数 中 加 入 邻 域 


别 、 医 学 影像 及 指纹 识别 等 各 种 图 像 处 理应 用 中 。 当 前 ， 


图 像 分 割 是 图 像 处 理 中 的 关键 所 在 ， 被 广泛 应 用 在 人 脸 识 项， 提出 了 基于 空间 约束 的 模糊 C- 均 值 算法 (FCM algorithm 
图 像 ”based on space constraint)， 增 强 了 抗 噪声 能 力 ， 但 也 增加 了 算 


分 割 方法 有 多 种 ， 比 如 区 域 增长 法 等 。 其 中 基于 聚 类 方法 被 广 ”法 的 时 间 复 杂 度 。 针 对 该 问题 ，Chen 等 人 6 用 邻 域 像素 的 均值 


mean 


司 , 人 硬 C- 均 值 聚 类 算法 以 及 模糊 C- 均 值 聚 类 算法 (fuzzyC- ”或 中 值 来 代 蔡 对 应 的 像素 点 提出 了 FCM_S1 和 FCM_S2 模型 


站 


s algorithm，FCMD)[I] 是 最 常用 的 两 种 聚 类 算法 。 与 HCM 该 模型 减少 了 算法 计算 时 间 , 但 模糊 了 图 像 的 许多 细节 。 于 是 ， 


算法 相 比 , FCM 算法 为 每 个 像素 的 归属 引入 了 模糊 概念 , 保留 。 Cai 等 人 0 提出 了 同时 考虑 邻 域 位 置 和 像素 关系 的 快速 FCM 


更 多 纹理 细节 信息 。 聚 类 算法 (fast generated fuzzy C-means algorithm )。 该 算法 虽 在 


相似 | 


生 ， 忽 视 了 图 像 其 他 影响 信息 ， 从 而 导致 FCM 的 分 割 结 ”不 高 。 
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对 于 参数 比较 多 同时 要 人 工 进行 设 定 问题 ，Stelios 等 个 
在 经 过 研究 之 后 ， 痔 述 了 局 部 空间 模糊 C- 均 值 算法 (fuzzy local 
information C-means algorithm，FLICM), 该 算法 引入 了 模糊 影响 
因子 ， 考 虑 邻 域 像 素 的 空间 信息 和 灰 度 信息 ， 且 参数 不 用 人 为 
设 定 , 大 大 增强 了 算法 适用 范围 ,其 全 部 参数 能 够 自 适应 确定 ， 
在 鲁 棒 性 方面 表现 非常 出 色 , 能 够 让 图 像 细 节 得 到 有 效 的 保存 。 
基于 此 ，Gong 等 人 在 局 部 空间 信息 模糊 C- 均 值 (fuzzy local 
information C-means algorithm) 算法 的 基础 上 引入 了 像素 点 的 
邻 域 像素 变异 系数 ck， 提出 了 修正 的 局 部 空间 模糊 C- 均 值 算 
法 (modified fuzzy local information C-means algorithm)。 若 噪 
声 高 度 污染 图 像 时 ， 在 一 个 像素 的 邻 域 像素 中 可 能 会 包含 很 多 
异常 特征 ， 如 图 1 所 示 ， 容 易 出 现 分 割 错 误 。 此 外 ，WELICM 
算法 用 欧 氏 距离 作 相 似 性 度量 ， 容 易 造 成 低 维 空间 线性 不 可 分 
现象 。 

为 解决 这 些 问 题 ， 本 文 在 模糊 影响 因子 中 引入 了 非 局 部 空 
间 信 息影 响 因 子 ， 根 据 图 像 噪声 污染 程度 同时 考虑 局 部 信息 和 
非 局 部 信息 对 中 心 像素 点 的 影响 ， 提 出 了 融合 局 部 和 非 局 部 空 
间 信 息 的 模糊 C- 均 值 聚 类 (NLWFLICMDJ) 算 法 ， 并 用 和 距离 代 蔡 
传统 欧 氏 距离 。 而 对 于 算法 本 身 依赖 于 初始 类 中 心 选择 和 复杂 
度 较 高 的 问题 ， 本 文 提出 了 Canopy-NLWFLICM 算法 ， 在 聚 类 
之 前 得 到 较为 准确 的 初始 类 中 心 , 在 此 基础 上 用 改进 算法 聚 类 ， 
实验 结果 也 说 明了 方法 的 有 效 性 及 对 实际 问题 的 实用 性 。 


1 ”相关 工作 


1.1 局 部 空间 信息 的 模糊 C- 均 值 聚 类 (WFLICM) 算法 

为 了 克服 传统 模糊 C- 均 值 聚 类 算法 在 图 像 分 割 时 仅 利 用 
像素 点 灰 度 值 信息 的 问题 ，WFLICM 算法 继承 了 FCM 算法 的 
优势 ， 同 时 还 融入 了 领域 空间 信息 ， 综 合 像素 点 邻 域 框 内 像素 
点 对 中 心 点 的 影响 ， 构 造 目标 函数 : 


SFalw ra] 0 


i=1 j= 


lk 


- 


其 中 :zi 在 WFLICM 算法 中 表示 第 个 像素 点 的 灰 度 值 , 戌 在 


WFLICM 算法 中 为 第 类 的 类 中 心 灰 度 值 ,至 j 是 加 入 的 模糊 影 


响 因子 ,相关 定义 具体 如 下 : 
1 2 
Em ,0 


其 中 : 本 表示 像素 点 区 所 在 的 邻 域 窗口 ; dd 表示 邻 域 窗口 中 
的 点 与 中 心 像素 点 的 坐标 位 置 距 离 ， 昭 ， 是 邻 域内 像素 
点 影响 因子 权重 ， 定 义 为 


1 
2+min((c./c) ,(c/ci) ) 
Wix = 1 = . 
,ife, <C 


2-minf(ejeJ Go 


,ife, >c 


陈 ” 凯 ， 等 : 基于 Canopy 聚 关 的 噪声 自 适应 模糊 C- 均 值 算 ; 


其 中 ; cy 是 邻 域 像素 点 的 变异 系数 ,以 中 心 像素 点 天 为 中 心 ， 
var (x) 
Hi 加 
其 中 ， Bar(x】 代表 邻 域 中 全 部 像素 点 灰 度 值 的 方差 , 下 是 邻 
域内 像素 点 灰 度 值 的 均值 。 


Cx 
E 二 (5) 


于 是 ， 由 利用 拉 格 朗 日 乘 子 法 得 到 隶属 度 更 新 公式 为 


1 . 
UR pl ENELEn 
ef d+G, 下 (6) 
大 =1 dx + Gx 
聚 类 中 心 更 新 公式 为 
Pus 
Wj= 芒 ,1<jsce (7) 
Su 


1.2 Canopy 算法 

Canopy 是 聚 类 算法 的 一 种 实现 , 它 是 一 种 快速 、 简 单 , 但 是 
不 太 精 确 的 聚 类 算法 , Canopy 聚 类 最 大 的 特点 是 不 需要 事先 指 
定 聚 类 的 个 数 ， 因 此 具有 很 大 的 实际 应 用 价值 ， 可 以 用 在 预 处 
理 阶 段 ， 对 数据 进行 粗 聚 类 。Canopy 算法 基本 流程 如 下 : 
a) 设置 阔 值 T1、T2, T1>T2 (关于 阔 值 的 选取 可 通过 交叉 
验证 得 到 ); 

b) 从 数据 集中 任 取 一 点 P， 将 P 作为 第 一 个 Canopy， 并 
将 它 从 数据 集中 删除 ; 

c) 继续 从 集合 中 取 点 ， 计 算 其 到 已 经 产生 的 所 有 Canopy 
的 距离 ， 如 果 到 某 个 Canopy 的 距离 小 于 T2， 则 将 其 加 入 到 该 
Canopy， 如 果 它 到 所 有 Canopy 中 心 的 距离 都 大 于 T1， 则 将 其 
作为 一 个 新 Canopy， 如 果 该 点 到 某 个 Canopy 距离 小 于 T1， 并 
在 其 与 所 有 Canopy 距离 计算 完成 后 仍 未 加 入 任何 Canopy， 则 
将 它 作 为 一 个 新 的 Canopy: 

d) 重复 第 c) 步 ， 直 至 数据 集 为 空 。 

算法 流程 图 如 图 1 所 示 。 


数据 集 Data， 疗 值 
Ta 


随机 选取 Data 中 P， 作 为 初始 
Canopy， 从 Data 中 删除 P 


继续 从 Data 中 取 点 ， 计 算 其 与 


已 有 canopy 间 的 距离 d 


vy 
将 所 有 d 大 于 T1 的 点 作为 新 
Canopy， 删 除 Data 中 小 于 T2 的 点 


图 1 Canopy 算法 流程 图 
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2 ”本 文 改进 算法 
2.1 融合 局 部 和 非 局 部 空间 信 
(NLWFLICM) 算 法 
WFLICM 算法 充分 考虑 了 局 部 空间 信息 和 局 部 灰 度 值 关 
系 对 像素 点 的 影响 ， 在 很 大 程度 上 降低 了 对 噪声 的 敏感 度 。 但 
当 图 像 噪声 污染 严重 时 ， 一 个 像素 点 的 邻 域 像素 会 存在 异常 特 
征 造成 局 部 信息 对 中 心 像素 点 聚 类 的 影响 发 生 错 误 。 对 于 这 种 
情况 ， 显 然 只 考虑 局 部 空间 信息 会 导致 图 像 分 割 效果 很 难 令 人 
满意 。 本文 提 出 融合 局 部 和 非 局 部 空间 信息 的 模糊 C- 均 值 聚 类 
算法 (NLWFLICM), 该 方法 将 非 局 部 空间 信息 加 入 到 影响 因 寺 


息 的 模糊 C- 均 值 聚 类 


¥O 
3 
je | 
X 
< 
np 
¥ 
es 
六 


陈 凯 ， 等 : 基于 Canopy 聚 类 


标准 差 , 它 可 以 通过 图 像 估算 得 到 [9], yr 是 图 像 的 平滑 参数 ， 
《本 文 F 设置 为 10)， 则 


刀 = ycr。 (13) 


模糊 因子 车; 中 控制 局 部 和 非 局 


部 空间 信息 的 权重 参数 @; 


可 用 定义 为 
4 = is) a 
其 中 Q 表示 搜索 窗口 内 所 有 像素 点 的 个 数 。 
由 于 使 用 相似 度 函 数 计算 权重 参数 ， 而 相似 度 函 数 取 决 于 


中 ， 并 根据 图 像 中 噪声 的 大 小 程度 在 去 噪 的 同时 尽量 保留 图 像 
的 细节 特征 。 定 义 目 标 函 数 (5) 中 的 模糊 影响 因子 为 


G; =(-wi 2 ee 四 ] -vf 


kzi 
keN; l (8) 


2 
+0, 5, (x,, 5h, ， 


rz#i 
reQ; 


1 
2 7 之 ee -wwP lr -wf, 
keN; 
Je， 加 »,) (10) 


> S;, (x,, 
rzi ‘ 
reQ, 


s 间 信息 影响 因子 ， 配 ;为 非 局 部 空间 信息 影响 


于 


则 向 为 局 部 空 


子 。 在 非 局 部 空间 信息 影响 因子 如; 中 ，S; (x, x 表示 搜索 


窗口 和 Q; 中 像素 x 与 中 心 像素 筷 间 的 相似 度 ,SG x) 满足 


条 件 0< Si(xi,xi) < 1 和 ZreniSir (Xr,Xi) = 1, 它 的 值 由 以 和.、 


Tr 

中 为 中 心 的 矩 邻 域 训 ， 与 下 间 的 距离 VCN) 一 VCNi) 上 决定 : 
vy -vy 

也 (11) 


ee 
“™ 3 


其 中 |IVCN,) 一 VCNi) 上 是 两 个 邻 域 中 像素 点 灰 度 值 均值 之 差 ， 


S(xi ) 为 轨 一 化 系数 : 


s(x,) De a ， (12) 


h 用 来 对 高 斯 函数 的 衰减 程度 进行 控制 。 h 的 大 小 影响 去 噪 效 
果 和 细节 的 保留 ， 越 大 去 品 能 力 会 更 强 ， 不 过 易 造 成 过 平滑 ， 
越 小 可 以 保留 更 多 图 像 细节 ， 不 过 去 噪 效果 会 变 差 。 本 文 参数 
可 以 根据 被 污染 图 像 的 噪声 级 别 自 适应 设置 。 设 本 是 噪声 的 


参数 h 的 值 ， 而 h 又 与 噪声 的 标准 差 成 正比 ， 所 以 当 噪 声 污染 
程度 较 小 时 ，h 的 值 将 会 很 小 ， 同 时 权重 参数 梧 也 会 很 小 ， 此 
时 主要 考虑 局 部 空间 信息 的 影响 ; 反之 ,， 当 噪声 污染 程度 较 大 
时 ， 非 局 部 空间 信息 对 中 心 像 素 点 聚 类 的 影响 较 大 。 因 此 根据 
自 适 应 参数 h ， 局 部 和 非 局 部 空间 信息 可 以 合理 地 得 到 利用 。 
本 文 用 高 斯 核 距离 ] 代 替 WFLICM 算法 中 传统 的 欧式 距离 。 
标 函 数 如 下 : 


l 


本 文 提出 NLWFLICM 算法 ，NLWFLICM 的 


G =(-w Em -us) le) - Bl) 
可 时 下 nt 


+Q; > 交配 让 


rzi 
reQ, 


2 3 oo 


了 


1 Gx) 一 (vw) 有 2 在 NLWFLICM 算法 中 表示 高 维特 征 空 


间 中 定义 的 距离 ，@() 则 表示 低 维 到 高 维 空间 的 上 映射。 具体 定 
义 为 
1 P00) — (wv) 吧 
= (eco — @(w)) (eco) 
一 ~(v)) 
= Px)T P(x) 一 (vw) Pxi) — P(x) TO(v;) 
+ P(v 站 ®(v;) 
K (xi, Xi) 中 K(v;, Dj) wm 2K(xi, Dj) 
(16) 
NLWFLICM 算法 中 核 函 数 使 用 的 是 高 斯 核 函 数 ; 
xly 
Klx,y) = cp 人 了 -) (17) 
于 是 ， 由 利用 拉 格 朗 日 乘 子 法 得 到 隶属 度 更 新 公式 为 
1 
Uij 二 N71) 
c (22x)+65 ) (18) 
Bs) 
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= Tii(uF kuvi)xi) 


(rk) U9 
2.2 Canopy-NLWFLICM 算法 
融合 局 部 和 非 局 部 空间 信息 的 模糊 C- 均 值 聚 类 算法 
(NLWFLICM) 解决 了 图 像 局 部 空间 信息 被 噪声 污染 导致 图 像 
分 割 精 度 下 降 和 分 割 错 误 的 问题 ， 但 在 算法 在 收敛 速度 和 算法 
复杂 度 上 和 传统 FCM 算法 一 样 依赖 初始 中 心 点 的 选取 ， 从 而 
收敛 速度 慢 等 缺点 。 为 此 本 文 在 


NLWFLICM 算法 聚 类 前 用 Canopy 算法 先 对 图 像 进 行 预 处 理 ， 


导致 算法 复杂 度 高 ， 


得 到 粗 聚 类 中 心 ， 提 出 Canopy-NLWFLICM 算法 ， 降 低 了 算法 
复杂 度 ， 提 高 了 聚 类 收敛 速度 。 

Canopy-NLWFLICM 算法 如 下 : 

a) 选 定 聚 类 数目 ce、 模糊 指数 mm、 邻 域 窗口 耻 、 搜 索 窗口 
、 最 大 迭代 次 数 T 和 闵 值 e、Tl 和 T2; 

b) 


Se 


| Canopy 算法 得 到 聚 类 初始 中 心 ; 


c) 利用 式 (4) 计 算 Wws ， 式 (1D 计 算 Si (xs 后) ， 


Sr 


d) 利用 式 (15) 计算 Gj; 


e) 利用 式 (18) 计 算 隶 属 度 和 矩阵 区 ， 利 用 式 (19) 计 算 聚 类 中 
心太 ; 

人) 根据 当前 矿 和 下 ， 计 算 目 标 函 数 (8) 的 值 。 若 迭 代 次 数 
大 于 T 或 者 相 邻 两 次 目标 函数 值 差 的 绝对 值 小 于 阔 值 在 ， 则 算 
法 停止 ， 否 则 ， 令 上 一 过 十 再 ， 转 b)。 


3 ”实验 结果 与 分 析 


为 了 测试 本 文 提出 的 Canopy-NLWFLICM 算法 图 像 分 割 性 
能 ,分别 在 人 工 合成 图 像 和 自然 图 像 上 进行 对 比 实验 。 
行 环境 是 Matlab 2014b。 以 下 实验 中 设 定 模糊 指数 m=2， 和 迭代 
闵 值 e=0.001， 最 大 迭代 次 数 T-500， 亚 为 3*3 邻 域 , 有; 为 
11*11 的 搜索 窗口 。 为 说 明 Canopy-NLWFLICM 算法 的 分 割 有 
效 性 ， 本 文 用 KFLICM、WFLICM 和 Canopy-FCM 作为 比较 算 


实验 运 


法 。 
3.1 人 工 合成 图 像 分 割 实验 

第 一 节 实 验 围绕 人 工 合成 图 像 进行 ， 总 数 为 3 张 ， 添 加 的 
高 斯 以 及 椒盐 噪声 强度 有 差异 。 对 比分 析 图 像 分 割 算 法 的 抗 品 
声 能 力 , 我 们 常 采用 分 割 准确 性 (SA) 指 标 来 衡量 其 好 坏 。SA 是 
指正 确 分 类 的 像素 点 占 总 分 割 像素 点 的 百分比 ， 代 表 图 像 分 割 


Chi 
4 


陈 ” 凯 ， 等 : 基于 Canopy 聚 关 的 噪声 自 适 应 模糊 C- 均 值 算 ; 


3 
兴 
a 
:> 
滞 
扎 一 
玉 


3 Ve, ， (19) 


j=l 


其 中 : c 为 分 割 的 类 数 ， 奋 ; 指 算法 分 割 到 第 类 的 像素 点 ， 
CC 指 图 像 中 属于 第 E 类 的 像素 点 全 

2 为 四 种 分 割 算 法 在 第 一 张 人 工 合成 图 像 中 添加 均值 为 
0 方差 为 0.02 高 斯 噪声 后 的 分 割 结果 。 从 分 割 结果 可 以 看 出 
KFLICM 算法 和 Canopy-FCM 算法 的 分 制图 受 噪 声 影响 非常 严 
重 ; WFLICM 算法 能 够 消除 一 定 的 噪声 ;， Canopy-NLWFLICM 
算法 在 去 噪声 和 保持 边缘 细节 上 表现 最 好 , 几乎 不 受 噪声 影响 ， 


表现 出 较 强 的 去 噪 能 力 。 
(c)KFLICM 分 制 区 


(9) 原 始 图 像 
| [A 


(DCanopy-NLWFLICM 

分 割 图 
图 2 第 一 张 合 成 图 像 分 割 实验 
基于 图 3 可 以 看 到 ， 随 着 高 斯 噪声 强度 越 来 越 大 ， 在 这 种 
情况 下 ， 对 比 算法 会 有 着 越 来 越 低 的 SA 值 ， 而 本 文 改 进 的 两 
种 算法 分 割 精 度 变化 最 小 ， 基 本 上 保持 在 99% 以 上 。SA 数据 
说 明 Canopy-NLWFLICM 算法 比 KFLICM、WFLICM 和 
Canopy-FCM 去 噪 效果 更 具 鲁 棒 性 。 


nt 
o 


(人 b) 噪 声 图 像 


(d)WFLICM 分 割 区 


(e)Canopy-FCM 分 割 医 
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图 3 四 种 算法 在 第 一 张 合成 图 像 上 的 分 割 精度 (SA%) 

第 二 张 合成 图 像 分 四 类 共 244*244 个 像素 点 ， 图 4 是 四 种 
算法 在 第 二 张 图像 受 密度 为 0.15 的 椒盐 噪声 污染 后 的 分 割 结 
果 。WFLICM 算法 受 噪声 影响 非常 严重 ， 分 割 图 上 留 有 大 量 的 
噪声 点 ; Canopy-NLWFLICM 算法 去 椒盐 噪声 能 力 最 强 ， 分 割 


结果 图 没有 明显 噪声 点 。 


As 
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(a) 原 始 区 


像 (人 b) 噪 声 


(d)WFLICM 分 割 


图 


图 4 第 二 张 合 成 


起 于 图 


情况 下 ， 对 比 算法 会 


种 算法 分 割 


第 三 张 合成 区 
算法 在 第 三 张 医 


显 下 降 ,，Canopy-NLWFLICM 
说 明 Canopy-NLWFLICM 
已 力 


噪声 能 


精度 变化 最 小 ， 其 


图 像 


(e)Canopy-FCM 分 割 


图 


分 割 


图 像 分 割 实验 


(c)KFLICM 分 割 区 


(f) Canopy-NLWFLICM 


区 


5 可 以 看 到 ， 随 着 椒盐 噪声 强度 越 来 越 大 ， 在 这 种 


有 着 越 来 越 低 的 SA 值 ， 而 本 文 改进 的 两 
他 两 种 种 算法 
算法 精度 下 降幅 度 最 小 。SA 数据 


的 分 割 精度 SA 明 


法 较 


其 他 三 种 算法 有 着 更 好 的 椒盐 


一 日 一 Canopy-FCM 
中 一 KFLICM 
92 吓 一 WFUICM 
一 一 Canopy-NLWFLICM 


(0.15) (020) 


图 5 


四 种 算法 在 第 - 


- 张 合成 


区 


像 分 四 


(025) 


(030) 


图 像 上 的 分 割 精度 (SA%) 


类 共 256*256 个 像素 点 。 


区 


6 为 四 


种 


区 


十 日 
结果 。 


NLWFLICM 算法 的 分 割 结 


KFLICM 


图 几乎 和 原 


像 添加 均值 为 0 方差 为 0.02 
算法 受 噪声 影响 严重 ， 分 割 医 
声 点 ; Canopy-NLWFLICM 算法 去 椒盐 噪声 和 


高 斯 噪声 后 的 分 割 


名 


上 留 有 大 量 的 噪 


图 


像 一 


E 力 最 强 
致 。 


RE, Canopy- 
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(a) 原 始 轿 


(b) 噪 声 区 


(CKFLICM 分 制图 


四 门 天 


(d)WFLICM 分 割 图 (e)Canopy-FCM 分 割 图 (DCanopy-NLWFLICM 
分 割 区 
图 6 第 三 张 合 成 图 像 分 割 实验 
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陈 凯 ， 等 : 基于 Canopy 聚 类 的 噪声 自 


基于 图 7 可 以 看 到 ， 随 着 高 斯 噪声 强度 越 来 越 大 ， 在 这 种 
情况 下 ， 对 比 算法 会 有 着 越 来 越 低 的 SA 值 ， 而 本 文 改 进 的 两 
种 算法 分 割 精度 变化 最 小 ， 其 他 三 种 算法 的 分 割 精度 SA 明显 
下 降 ， Canopy-NLWFLICM 算法 精度 下 降幅 度 最 小 。SA 数据 
说 明 本 文 改进 的 算法 较 其 他 三 种 算法 去 噪 能 力 更 强 。 


图 7 四 
3.2 自然 图 像 分 割 实验 

本 实验 用 三 张 自 然 图 像 来 比较 
图 像 时 的 分 割 性 能 。 在 比较 四 种 算法 
力 时 将 涉及 到 结构 相似 性 (SSIM),SSIM 
似 性 程度 : 


种 算法 在 第 二 


图 像 上 的 分 割 精度 (SA%o) 


0 种 算法 的 在 处 理 自然 噪声 
在 自然 图 像 上 的 
j 来 判断 两 幅 


Qn, + clocc， 十 cja 十 oj 


SSIM = 对 村 六 + ce: 证 0 + cjec 小 Cj) 


(20) 


其 中 


C=(KD’,C, =(K,D)’,Cs = 字 ,K =0.01 ， 


K, =0.03,L=255 ， 平 均值 和 方差 的 计算 公式 如 下 : 


CD) 


m An 


4), (23) 


Oxy = 


EEC- Je 让- 

isl 1 
此 外 实验 采用 欠 分 割 率 CUR) 对 分 割 

UR = 二 

R. + 0O, 
其 中 : 栈 , 表示 本 应 该 包含 在 分 割 结果 中 实际 却 不 在 分 割 结果 中 
的 像素 点 个 数 ， 民 ,表示 理论 分 割 图 像 中 像素 点 个 数 ， 人 表示 

不 应 该 包含 在 分 割 图 像 中 的 像素 点 个 数 。 

第 一 张 自 然 图 像 大 小 分 别 为 173*173， 图 8 给 出 了 四 种 算 
法 在 第 二 张 图 像 添加 椒盐 (0.10) 品 声 后 的 分 割 结果 。 从 第 一 张 
然 图 像 分 割 效果 图 明显 可 以 看 出 , WFLICM 算法 出 现 了 分 割 错 
误 ， Canopy-FCM 算法 和 KFLICM 算法 边缘 分 割 不 够 清晰 , 噪 
明显 ， 而 Canopy-NLWFLICM 算法 分 割 图 最 接近 原 网 。 


结果 进行 评价 。 


(24) 


201804.02368v1 


chinaXiv 
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(c)KFLICM 分 割 


(b) 噪 声色 


(d)WFLICM 分 割 图 (e) Canopy-FCM 分 割 图 


(fCanopy-NLWFLICM 
分 割 图 
图 8 第 一 张 自 然 图 像 分 割 实验 


表 1 给 出 了 四 种 算法 分 割 结果 图 与 第 一 张 自 然 图 像 间 的 结 
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构 相 似 性 。SSIM 的 数据 说 明 Canopy-NLWFLICM 分 制图 与 原 
图 结构 上 最 为 相似 。 
表 3 四 种 算法 分 割 图 与 第 二 张 自 然 图 像 间 的 结构 相似 性 


噪声 KFLICM WFLICM CanopyFCM Canopy-NLWFLICM 
椒盐 (0.10) 0.9335 0.9422 0.9366 0.9597 
椒盐 (0.15) 0.9283 0.9318 0.9295 0.9536 
椒盐 (0.20) 0.9107 ”0.9259 0.9131 0.9498 


表 4 给 出 了 四 种 算法 在 第 一 3 
据说 明 Canopy-NLWFLICM 算法 在 分 割 ; 


然 图 像 上 的 欠 分 割 率 ， 数 
佳 确 率 相 比 另 三 种 算法 


表 4 四 种 算法 在 第 二 张 自 然 图 像 上 的 欠 分 割 率 (UR%) 


加 噪声 。 KFLICM WFLICM CanopyFCM Canopy-NLWFLICM 
DN DBT SE Canopy ETOM oo ad ea 0.0919 
图 结构 上 最 相似 。 椒盐 (0.20) 0.1977 0.1536 0.1953 0.1232 
表 1 四 种 算法 分 割 图 与 第 一 张 日 然 图 像 间 的 结构 相似 性 椒盐 (0.25) 0.2353 0.1877 0.2322 0.1604 
噪声 。 KFLICM WFLICM CanopyFCM Canopy-NLWFLICM 
椒盐 (0.10) 0.9515 0.9428 0.9503 0.9769 elon ni 
法 在 第 二 张 图 像 添加 椒盐 (0.02) 品 声 后 的 分 割 结果 。Canopy- 
0 a NLWFLICM 分 割 效果 图 比 其 他 三 种 算法 分 割 效果 图 去 噪 效果 
椒盐 (0.20) 0.9301 0.9243 0.9274 0.9682 


表 2 给 出 了 四 种 算法 在 第 一 张 自 然 图 像 上 的 从 分 割 率 ， 数 
据说 明 Canopy-NLWFLICM 算法 在 分 割 准确 率 相 比 另 三 种 算法 


表 2 四 种 算法 在 第 一 张 自 然 图 像 上 的 欠 分 割 率 (UR%) 


噪声 。 KFLICM WFLICM CanopyFCM Canopy-NLWFLICM 
概 盐 (0.10) ”0.1856 0.1271 0.1583 0.0774 
椒盐 (0.15) ”0.2539 0.1694 0.1843 0.1109 
椒盐 (0.20) ”0.3328 0.2208 0.2535 0.1592 


第 二 张 自 然 图 像 大 小 分 别 为 173*173， 图 9 给 出 了 四 种 算 
法 在 第 二 张 图 像 添加 椒盐 (0.15) 噪 声 后 的 分 割 结果 。Canopy- 
然 图 像 上 去 椒盐 噪声 能 力 最 差 ， 
NLWFLICM 分 割 效果 图 比 另 三 种 算法 在 去 噪 和 保留 图 像 细 节 
上 的 效果 更 好 。 


心 


FCM 算法 在 


Canopy- 


,外 原始 图 像 


(d)WFLICM 分 割 图 (e) Canopy-FCM 分 割 图 


(fjCanopy-NLWFLICM 
分 割 图 
图 9 第 二 张 自 然 图 像 分 割 实验 


表 3 给 出 了 四 种 算法 分 割 结果 图 与 第 二 张 自 然 图 像 间 的 结 


更 好 ， 图 像 背景 与 人 物 细节 分 割 最 为 清晰 ， 也 更 加 光滑 流畅 。 


(e)Canopy-FCM 分 割 图 


(d)WFLICM 分 割 图 


(fCanopy-NLWFLICM 
分 割 区 
图 10 第 三 张 自 然 图 像 分 割 实 验 


表 5 给 出 了 四 种 算法 分 割 结果 图 与 第 三 张 自 然 图 像 间 的 结 
构 相 似 性 。SSIM 的 数据 明显 说 明 Canopy-NLWFLICM 分 割 
与 原 图 结构 性 最 相似 。 

表 5 四 种 算法 分 割 图 与 第 三 张 自 然 图 像 间 的 结构 相似 性 


噪声 KFLICM WFLICM CanopyFCM Canopy-NLWFLICM 
椒盐 (0.01) 0.8925 0.8742 0.8906 0.9039 
椒盐 (0.02) 0.8276 ”0.7906 0.8538 0.8577 
椒盐 (0.03) ”0.7503 0.7255 0.8322 0.8407 


表 6 给 出 了 四 种 算法 在 第 一 张 自 然 图 像 上 的 欠 分 割 率 ， 数 
据说 明 Canopy-NLWFLICM 算法 在 分 割 准确 率 相 比 男 三 种 算法 


表 6 四 种 算法 在 第 三 张 自 然 图 像 上 的 欠 分 割 率 (UR%) 


噪声 KFLICM WFLICM CanopyFCM Canopy-NLWFLICM 
椒盐 (0.10) 0.1905 0.1866 0.1952 0.1644 


录用 稿 
椒盐 (0.15) ”0.2283 ”0.2103 0.2313 0.1893 
椒盐 (0.20) 0.2579 ”0.2399 0.2679 0.2271 


以 上 所 有 试验 都 证 明 Canopy-NLWFLICM 算法 比 KFLICM 
算法 和 WFLICM 算法 有 更 强 的 去 噪声 能 力 , 相 比 Canopy-FCM 
算法 能 更 好 地 保存 图 像 细 节 ， 与 原 图 的 结构 性 也 最 为 相似 。 此 
外 ， 在 大 量 噪声 污染 时 ，WFLICM 算法 出 现 分 割 错误 问题 ， 而 
Canopy-NLWFLICM 算法 会 通过 噪声 程度 对 局 部 和 非 局 部 信息 
进行 对 应 地 调整 ， 从 而 克服 分 割 错误 所 带 来 的 问题 。 


4 ”结束 语 


针对 局 部 空间 信息 的 模糊 C- 均 值 算法 (WFLICM) 在 被 大 
量 噪声 污染 时 分 割 错 误 问 题 ， 本 文 引 入 非 局 部 信息 影响 因子 ， 
提出 了 NLWFLICM 算法 , 实验 结果 表明 该 方法 具有 更 好 的 图 
力 。 由 于 参数 h 可 以 根据 图 像 中 存在 的 噪声 标准 差 计 
使 其 适应 于 噪声 的 比例 ， 而 局 部 与 非 局 部 空间 的 权重 可 通 
过 非 局 部 空间 信息 来 计算 ， 从 而 使 局 部 和 非 局 部 空间 信息 自 适 
应 地 影响 目标 函数 中 的 相似 性 度量 。 以 这 种 方式 ， 局 部 和 非 局 
部 空间 信息 可 以 被 适当 地 用 于 低 程 度 和 高 程度 损坏 的 噪声 图 像 
的 分 割 。 针 对 算法 依赖 初始 类 中 心 选择 问题 , 本 文 提出 Canopy- 
NLWFLICM 算法 ， 首 先 用 Canopy 算法 进行 预 处 理 ， 得 到 初始 
类 中 心 ， 在 此 基础 上 用 NLWFLICM 算法 进行 聚 类 ， 提 高 了 算 
法 的 分 割 精度 和 收敛 速度 。 通 过 不 同类 型 的 分 割 实 验证 明了 本 
文 算法 的 优秀 的 去 噪声 的 能 力 和 良好 的 鲁 棒 性 ， 图 像 细节 也 同 
时 得 到 一 定 的 保留 ， 具 有 一 定 的 实际 应 用 价值 。 当 然 ， 本 文 也 
存在 一 些 不 足 之 处 ， 其 算法 复杂 度 较 高 ， 接 下 来 的 工作 就 是 在 
此 基础 上 提高 时 间 效 率 ， 使 其 更 具 实 用 性 
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